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ABSTRACT 

Upcoming 7-ray satellites will search for Dark Matter annihilations in Milky Way 
substructures (or 'clumps'). The prospects for detecting these objects strongly depend 
on the assumptions made on the distribution of Dark Matter in substructures, and on 
the distribution of substructures in the Milky Way halo. By adopting simplified, yet 
rather extreme, prescriptions for these quantities, we compute the number of sources 
that can be detected with upcoming experiments such as GLAST, and show that, 
for the most optimistic particle physics setup (m x = 40 GeV and annihilation cross 
section av = 3 x 10~ 26 cm 3 s _1 ), the result ranges from zero to ~ hundred sources, all 
with mass above 1O 5 M0. However, for a fiducial DM candidate with mass m x = 100 
GeV and av — 10~ 26 cm 3 s -1 , at most a handful of large mass substructures can be 
detected at 5a, with a 1-year exposure time, by a GLAST-like experiment. Scenarios 
where micro-clumps (i.e. clumps with mass as small as 10 _6 M©) can be detected are 
severely constrained by the diffuse 7-ray background detected by EGRET. 



1 INTRODUCTION 

Indirect Dark Matter [DM] searches are based on the detec- 
tion of secondary particles and radiation produced by the 
self-annihilation of DM particles (Bergstrom 2000, Bertone 
et al. 2005a). 

Although the predicted annihilation flux is typically af- 
fected by large astrophysical uncertainties, the detection of 
multiwavelenght photons, neutrinos or anti-matter from re- 
gions with high DM density would be of paramount impor- 
tance for the identification of DM particles. In fact, accel- 
erator searches of Physics beyond the Standard Model at 
the Large Hadron Collider, will not necessarily unveil the 
nature of DM, even if new particles are discovered, due to 
the difficulties associated with the reconstruction of the cos- 
mological abundance of the newly discovered particles (e.g. 
Baltz et al. 2006a, Nojiri et al. 2005). At the same time, 
DM particles could have small enough couplings to nucle- 
ons, to lead to null searches in direct detection experiments 
(see e.g. Munoz 2003 and references therein). 

In the framework of indirect DM searches, several 
strategies have thus been devised, in order to obtain conclu- 
sive evidence from astrophysical observations. For instance, 
one could search for peculiar features, such as lines or sharp 
cut-offs, in the 7-ray spectrum. Although for commonly 
studied DM candidates there are no tree level processes for 
direct annihilation into photons, loop-level processes to 77 
and 7Z may produce detectable lines at an energy equal to 



the DM particle mass (see e.g. Bergstrom & Ullio 1997, Ullio 
& Bergstrom 1998, Gounaris et al. 2003, Bergstrom et al. 
2005a). Other spectral features may help distinguishing the 
DM annihilation signal from ordinary astrophysical sources 
(Bergstrom et al. 2005b, Bergstrom et al. 2005c; see also 
the discussion in Baltz et al. 2006b). Alternatively, one can 
search for annihilation radiation from regions characterized 
by large concentrations of DM, but very few baryons, such 
as DM substructures in the Milky Way [MW hereafter] halo, 
including dwarf galaxies (Baltz et al. 2000, Tasitsiomi et 
al. 2003, Pieri & Branchini 2003, Evans et al. 2003, Tyler 
2002, Colafrancesco et al. 2007, Bergstrom & Hooper 2006) 
and DM mini-spikes around Intermediate Mass Black Holes 
(Bertone 2006, Bertone et al. 2005b, Horiuchi & Ando 2006, 
Fornasa et al. 2007, Brun et al. 2007). Finally, DM an- 
nihilation features can be detected in the energy spectrum 
and angular distribution of the cosmic 7-ray background 
(Bergstrom, Edsjo & Ullio 2001, Ando & Komatsu 2006). 

In the popular Cold DM scenario, gravitational insta- 
bilities lead to the formation of a wealth of virialized struc- 
tures, the DM haloes, spanning a huge range of masses, from 
the largest clusters of galaxies of ~ 1O 15 M0 down to Earth- 
size clumps of ~ 1O _6 M0 (Green et al. 2004, Green et 
al. 2005). Although the detectability of individual DM sub- 
structures, or "clumps" has been widely discussed in liter- 
ature, the number of detectable clumps with a GLAST-like 
experiment, at 5<r in 1 year and for a WIMP DM parti- 
cle is highly uncertain, ranging from <, 1 (Koushiappas 
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et al. 2004) to more than 50 (Baltz 2006b) for large mass 
haloes, while for microhaloes (i.e. clumps with a mass as 
small as 10 _6 M©) the predictions range from no detectable 
objects (Pieri et al. 2005) to a large number of detectable 
objects, with a fraction of them exhibiting a large proper 
motion (Koushiappas 2006). The apparent inconsistency of 
the results published so far, is actually due to the differ- 
ent assumptions that different groups adopt for the physical 
quantities that regulate the number and the annihilation 
"brightness" of DM clumps. In particular, even in the con- 
text of the benchmark density profile introduced by Navarro, 
Frenk and White 1996 [NFW], the results crucially depend 
on the substructures mass function, their distribution within 
the halo host and their virial concentration c(M, z) which is 
a function of mass and of collapse redshift of DM clumps. 

The paper is organized as follows: in Sec. 2, we describe 
the model we have adopted for the smooth component of 
the Galactic halo, and introduce the eight different models 
for DM substructures that will be discussed in the rest of 
the paper. In Sec. 3 we estimate the contribution to the 7- 
ray flux due to the smooth Galactic halo, unresolved DM 
clumps, and resolved (detectable) clumps. In Sec. 4, we study 
the prospects for detection of substructures with upcoming 
experiments such as GLAST, and in Sec. 5 we discuss the 
results and present our conclusions. 



2 MODELING THE GALACTIC HALO AND 
ITS SUBSTRUCTURES 

High resolution N-body experiments indicate that a large 
fraction of the mass within Dark Matter haloes is in the form 
of virialized subhaloes in all resolved mass scales. Their an- 
nihilation signal, which adds to that of the smooth Galactic 
component could be significant (Stoher et al. 2003, Die- 
mand et al. 2006, Diemand et al. 2007a). A precise model- 
ing of both the smooth DM distribution (the diffuse galactic 
component) and the subhalo population within (the clumpy 
component) is therefore mandatory to assess the possibility 
of detecting DM subhaloes through their annihilation sig- 
nal. High resolution numerical simulations enable to study 
gravitationally bound subhaloes with Msh > 10~ 6 MHaio, 
where Afnaio is the mass of the host, and therefore cannot 
resolve substructures in a MW-size halo all the way down 
to ~ 10~ 6 M©. In fact, the smallest substructures of inter- 
est must be studied within host haloes with mass ~ O.IMq 
and only at very large redshift (Diemand et al. 2006). At 
z — and within a MW-host the smallest subhaloes that 
we resolve have masses > 10 6 M Q (Diemand et al. 2007a) 
As a consequence the spatial distribution, mass function 
and internal structure of Galactic subhaloes can only be 
interpolated from the results of several numerical experi- 
ments spanning a large range of masses and redshifts using 
self-similarity arguments. This interpolation is affected by 
a number of uncertainties that we account for by exploring 
different models that meet the numerical constraints. 

2.1 The diffuse Galactic component 

The recent "Via Lactea" high resolution simulation (Die- 
mand et al. 2007a) shows that the density profile of a MW- 
sized DM halo is consistent to within 10% with the NFW 



profile that we adopt here: 



(1) 

where r is the distance from the halo center. This profile 
depends on two free parameters, the scale density, p s , and 
the scale radius, r 3 , that are related to each other by the 
virial mass of the halo, Mh ■ The latter is the mass enclosed 
in a sphere with radius r vir within which the mean density is 
200 times above critical. A different definition of r„; r would 
not change the DM profile nor our results. . Finally, we adopt 
a small core radius r min — 10 kpc. 

An important shape parameter that characterize the 
density profile is the virial concentration defined as the 
ratio between the virial radius and the scale radius, c = 
r v i r (Mh)/r s . Theoretical considerations corroborated by nu- 
merical experiments show that a relation exists between the 
mass of a halo, its collapse redshift, z co ii, and concentration 
parameter c. The collapse redshift is defined as in Bullock et 
al. 2001, as the epoch in which a mass scale Mh breaks into 
the nonlinear regime, i.e. when a(Mh)D(z co ii) ~ 1, where 
a(M h ) is the present linear theory amplitude of mass fluc- 
tuation on the scale Mh and D(z co n) is the linear theory 
growth factor at the redshift z co \\. The two models proposed 
by Bullock et al. 2001 [B01] and Eke et al. 2001 [ENS01] 
are consistent for masses larger than ~ 1O 9 M0 and in this 
paper we use the concentration parameter by B01 to model 
the diffuse Galactic component. On the contrary, for smaller 
masses ENS and B01 predictions become very different and 
we will have to consider both of them to model the subhaloes 
annihilation signal. 

Since a sizable fraction of the mass in the MW is in 
form of virialized subhaloes, there is not a unique way to 
determine p s , and r s of the host MW halo. For this rea- 
son we adopted two different procedures. In the first one 
we have computed p s , and r s as if the total mass of the sys- 
tem, including the clumpy component, were in fact smoothly 
distributed in the Galactic halo. In the second we have con- 
sidered the total mass of the system to determine the con- 
centration parameter but have used the diffuse component 
alone to relate p s to r s . Having checked that both proce- 
dures give similar predictions for the probability of subhalo 
detection (Section 3), in the following we will only discuss 
models based on the second approach. For the Milky Way we 
have used a virial mass Mh = 10 12 M Q and a concentration 
parameter c„; r ~ 9.8. 



2.2 The clumpy component 

To account for the presence of a population of DM subhaloes 
and investigate their effect on the annihilation signal, we 
need to specify their mass spectrum, spatial distribution and 
density profile. 
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Table 1. Halo model parameters. Column 1: Model name. Col- 
umn 2: Halo density profile for M > 1O 4 M0. Column 3. Halo 
density profile benchmark for M = 1O -6 M0, used as normaliza- 
tion. Column 4: Collapse rcdshift model, used as normalization. 
Column 5: Concentration parameter at 10 — ®Mq. For all the mod- 
els we have adopted an NFW profile. 

2.2.1 The Mass function and the Spatial Distribution of 
Clumps 

High resolution N-body experiments show that the mass 
function of both isolated field haloes and subhaloes is well 
approximated by a power law 

dn(M)/dln(Af) oc M~ a , (2) 

with a = 1, independently of the host halo mass, over 
the large redshift range z = [0,75], and mass intervals 
M = [Hr 6 ,lO lo ]M (Jenkins et al. 2001, Moore et al. 
2001 Diemand et al. 2004, Gao et al. 2005, Reed et al. 
2005, Diemand et al. 2006). Self-similarity is preserved at 
the present epoch down to the smallest masses if subhaloes 
survive gravitational disturbances during early merger pro- 
cesses and late tidal disruption from stellar encounters. 

Analytical arguments have been given against (Zhao et 
al. 2005) or in support of this hypothesis (Moore et al. 2005, 
Brezinsky et al. 2006). Moreover, currently resolved mass 
functions in numerical experiments suffer from dynamical 
friction at the high mass end which could steepen the halo 
mass function. Since changing the halo mass function slope 
might have a non-negligible impact on our analysis, we adopt 
a power-law index a — 1 as a reference case but also explore 
two shallower subhalo mass functions with a — 0.9 and a = 
0.95. All plots in this paper refer to the reference case a = 1 
and discuss the effect of adopting shallower slopes in the 
text. 

As far as the spatial distribution of subhalos inside our 
Galaxy is concerned, we follow the indications of the numer- 
ical experiment of Reed et al. 2005 and assume that the 
subhalo distribution traces that of the underlying host mass 
from r v ir and down to a minimum radius, r min (M), within 
which subhaloes are efficiently destroyed by gravitational 
interactions. We explicitly assume spherical symmetry and 
we ignore the possibility, indicated by some numerical ex- 
periments, that the radial distribution of subhalos might be 
more extended than that of the dark matter. 



Folding these indications together we model the number 
density of subhaloes per unit mass at a distance R from the 
GC as: 

Psh (M, R) = AM~> {R /J^]^^L )2 Me 1 kpc' 3 ,(3) 

where r^' IW is the scale radius of our Galaxy and the ef- 
fect of tidal disruption is accounted for by the Heaviside 
step function 6(r — r min (M)). To determine the tidal ra- 
dius, r m in{M), we follow the Roche criterion and compute 
it as the minimum distance at which the subhalo self gravity 
at r s equals the gravity pull of the halo host computed at 
the orbital radius of the subhalo. As a result r min (M) is an 
increasing function of the subhalo mass, implying that no 
subhaloes survive within r m i„(10 -6 Mq) ~ 200 pc. 

To normalize eq. 3, we again refer to numerical simula- 
tions that show that [5-10]% of the MW mass is distributed 
in subhaloes with masses in the range 10 7 - 1O 1O M (Die- 
mand et al. 2005, Diemand et al. 2007) In the following 
we use the optimistic value of 10 % and note that assuming 
5%, instead, would decrease the probability of subhalo de- 
tection by a factor 2. With this normalization about 53% of 
the MW mass is condensed within ~ 1.5 x 10 16 subhaloes 
with masses in the range [10~ 6 , 1O 1O ]M0, whose abundance 
in the solar neighborhood is remarkably high (~ 100 pc -3 ). 
The remaining 47% constitutes the diffuse galactic compo- 
nent that is assumed to follow a smooth NFW profile. We 
do not account here for the presence of mass debris streams 
resulting from tidal stripping since these structures, char- 
acterized by a mild density contrast, would not contribute 
significantly to the annihilation flux. 



2.2.2 Density Profile. 

Finally, we need to specify the density profile for the sub- 
structures. Constraints from numerical models only applies 
to limited mass ranges at very different epochs. At z = Die- 
mand et al. 2007b find that the velocity profile of Galactic 
subhaloes above 4 x 10 6 Mq in the "Via Lactea" simulation 
are well fitted by the NFW model. This result is also valid 
for the much smaller substructures with masses in the range 
[10~ 6 , 4 x 10~ 3 ] M Q that populate a parent halo of 0.014M© 
at 2 = 86 (Diemand et al. 2006). A large fraction of these 
small substructures do not survive the early stage of hierar- 
chical merging and late tidal interaction with stellar encoun- 
ters (Zhao et al. 2005). The ~ 10 16 survivors suffer from 
significant mass loss. Presumably this modifies their original 
NFW density profile that, however, seems to be preserved in 
the innermost region where most of the annihilation signal 
originates from (Kazantzidis et al. 2004). 

Yet these constraints from numerical experiments do 
not uniquely define the subhalo density profiles. Therefore, 
instead of relying on a single model profile we will explore 
several of them in an attempt of bracketing the theoretical 
uncertainties. 

All models that we have considered, and that are listed 
in Table 1, assume that subhaloes have the same NFW den- 
sity profile as their massive host but with different concen- 
tration parameters. These models have been flagged with 
the following prescriptions: 
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• B- models assume the c(M) relation of B01 for M > 
1O 4 M ; 

• SA^-models, use the ENS01 model for M > 10 4 M Q ; 

• In models flagged as z$ the low mass extrapolation of 
the concentration parameter is normalized to that of field 
haloes of 10 -6 M© measured in the numerical simulation of 
Diemand et al. 2005 (DMS05) at z = 26 and linearly ex- 
trapolated at z = 0, i.e. c(10 _6 M Q ,z = 0) = c(lO _6 M ,z = 
26) x (1 + 26). The underlying assumption is that, as in 
the Press- Schechter approach, all existing haloes have just 
formed, and the (1 + 26) scaling is required to account for 
the change in the mean density between z = 26 and z = 0; 

• Models flagged as z c assume instead that, once formed, 
subhaloes that survive to z = do not change their density 
profile. Therefore for each subhalo of mass M we determine 
its collapse redshift z c (defined as in B01 for M > 1O 4 M ) 
and compute its concentration parameter accordingly, as 
c(M,z = z c ) = c(M,z = 0)/(l + z c ), where c(M,z = 0) 
is computed as in the zO case. The small mass normaliza- 
tion is the same as in models zO, rescaled for the smallest 
masses collapse redshift z c (M — 1O _6 M ) — 70 suggested 
by DMS05 and used by Koushiappas 2006. In Fig. 1 we show 
the collapse redshift as a function of the mass adopted in the 
z c -models (thin line); 

• Models flagged 5a assume that all existing subhaloes 
with mass M = 1O _6 M form at the 5a peaks of the density 
field, which has the effect of increasing the normalization to 
c(10~ 6 M o ,z = 0) = 400; 

• Models flagged as ref, zO and ref, zc are identical to 
the corresponding zO and zc models for M > 1O 4 M0. In- 
deed, they use a naive extrapolation of the B01 model at 
low masses, both for the collapse redshift (thick line in Fig. 

1) and for the concentration parameter (filled circles in Fig. 

2) . Though this extrapolation is not supported by numerical 
simulations, it intuitively reflects the theoretcal flattering of 
the a(M) curve at low masses. 

The c(M) profile for each of the zO models is shown in Fig. 2. 

Finally, as pointed out in B01, the c(M) relation is not 
deterministic. Instead, for a fixed mass, the probability of a 
given value for c(M) is well described by a lognormal distri- 
bution 



P(c(M)) 



2ira c c(M) 



(c(M))-ln(c(M))' 



(4) 



where the mean c = c(M) is the concentration parameter of 
B01 or ENS and the dispersion a c = 0.24 does not depend 
on the halo mass (B01). In the following, we include this 
lognormal scatter in all models described in Table 1. 

We will not consider in this paper the possibility that 
subhaloes might contain sub-substructres (Strigari et al. 
2007, Diemand 2007a) and that their concentration (and 
scatter) might depend on the distance from the Galactic 
Centre (Diemand et al. 2007b). Indeed, such features are 
found in numerical simulations capable of resolving large 
sub-haloes (M > 1O 6 M ) but there is no evidence whether 
they also apply to much smaller haloes, with masses down 
to 1O~ 6 M . 



100 




Halo Mass [M„ n ] 

Figure 1. Collapse redshift, 2 C , as a function of the halo mass 
for model B01 (points) for M > 1O 4 M . The thin line shows the 
extrapolation at low masses normalized to the DMS05 value for 
2 C (10 — ®Mq). The thick line shown a naive extrapolation of the 
B01 model at low masses. 



3 7-RAY FLUX FROM ANNIHILATION IN DM 
CLUMPS 

The photon flux from neutralino annihilation in galactic sub- 
haloes, from a direction in the sky making an angle ip from 
the Galactic Center (GC), and observed by a detector with 
angular resolution 6, can be factorized into a term depending 
only on particle physics parameters, d$ pp / dE 1 and a term 
depending only on cosmological quantities, $ cosmo (-!/>, #): 



IK 



(£ 7 ,V,0) 



dE~ 



-(E-y) x $ cosmo (V,6>) 



(5) 



3.1 Particle physics contribution 

The first factor of Eq.5 can be written as: 



d& 



dE~ 
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47T 2m| 



(6) 



where m x denotes the Dark Matter particle mass and 
dN^/dE^ is the differential photon spectrum per annihi- 
lation relative to the final state /, with branching ratio 
Bf. Although the nature of the DM particle is unknown, 
we can make an educated guess on the physical parame- 
ters entering in the above equation. The most commonly 
discussed DM candidates are the so-called neutralinos, aris- 
ing in Super-symmetric extensions of the Standard Model 
of particle physics [SUSY], and the B^ particles, first ex- 
citation of the hypercharge gauge boson in theories with 
Universal Extra Dimensions [UED] (see Bergstrom 2000, 
Bertone et al. 2005b and references therein). Typical val- 
ues for the mass of these candidates range from ~ 50 GeV 
up to several TeV. The annihilation cross section can be as 
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Figure 2. Concentration parameters as a function of halo 
mass at z = computed for the ENS Z0 (solid), B zo (dotted) 
and the B Z „5 CT (dashed) model described in the text. Filled cir- 
cles show the B01 model naively extrapolated at low masses 
(B re/:Z0 model). 



Figure 4. $ cosmo a s a function of the angular distance from 
the Galactic center, for the MW smooth component (solid thick 
line) and for unresolved clumps in model B Zej 5 CT (solid thin line), 
B Zc (dot-dashed), ENS Zc (long dot-dashed), B ZOj 5 CT (long dashed), 
B re ^ Zc (open circles), B zo (dotted), B re ^ zo (filled circles) and 
ENS ZQ (dashed). The mass function slope is a = 1. 
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Figure 3. Exclusion plot in the (u anI1 !), «^x) pl ane - The solid line 
corresponds to our best case particle physics scenario, adopted 
for all the zO models. Models above the dashed (long dot-dashed, 
dot-dashed) curve violate the EGRET extra galactic background 
(EGB) flux constraints in scenario ENS Zc (B Zc ,B Zct 5 rT ) for a mass 
function slope a = 1 . For comparison, we show with filled (empty) 
dots a scan of SUSY models with relic density within 2 (5) stan- 
dard deviations from the WMAP+SDSS value. The dotted line 
corresponds to the UED model. Supersymmetric models were ob- 
tained using the DARKSUSY package (Gondolo et al. , 2004) 



high as CTannf = 3 x lCP 26 cm a s~ 1 , as appropriate for ther- 
mal relics that satisfy the cosmological constraints on the 
present abundance of Dark Matter in the Universe. How- 
ever, we note that the annihilation cross section can be much 
smaller, as the appropriate relic density can be achieved 
through processes such as co- annihilations (Bergstrom 2000, 
Bertone et al. 2005a). This can be seen from Fig. 3 that 
shows the range of o" ann u and m x allowed in the UED and 
SUSY models: solid (empty) circles correspond to models 
with relic density within 2 (5) standard deviations from 
the WMAP+SDSS suggested value Q, DM h 2 = OTOSOl^J, 
where as usual Q.dm is the DM density in units of the crit- 
ical density, and h is the Hubble parameter in units of 100 
km s- 1 kpc- 1 , h = 0.730±E!'^ (Tegmark et al. 2006). 

In order to optimize the prospects for detection, we 
adopt here a very low value for the particle mass, m x = 40 
GeV, together with a high annihilation cross section <r a nn« = 
3 x 10 _26 cm 3 s _1 . As for the nature of the DM particle, we 
assume here a 100% branching ratio in bb and the dN^ b /dE 7 
functional form of Fornengo et al. 2004. The results can be 
rescaled for any other candidate, although in most cases the 
photon spectrum arising from annihilations yields similar 
results. 



3.2 Annihilation flux from diffuse matter and 
unresolved clumps 

The contribution of unresolved substructures to the annihi- 
lation signal is given by 

3> cosmo (V>,Afi) = [ dM I dc ( ( d6d<j) ( d\ 

J M Jc J J A!Tl Jl.o.s 
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jcosmo 



[p sh {M,R(Re,\,iP,e,cj>)) x P(c)x 
(M,c,r(X,X',^,e', <!>')) x J(x,y,z\X,6,4>)] 



(7) 



where Af2 is the solid angle of observation pointing in the 
direction of observation ip and defined by the angular res- 
olution of the detector 0; J(x,y, z\X,9, (p) is the Jacobian 
determinant; R is the galactocentric distance, which, inside 
the cone, can be written as a function of the line of sight (A) 
and the solid angle (6 and <j>) coordinates and the pointing 
angle ip through the relation R — \J X 2 + Rq 2 — 2XRqC, 
where Rq is the distance of the Sun from the Galactic Cen- 
ter and C — cos(#) cos(V') — cos(0) sin(#) sin(-0); r is the 
radial coordinate inside the single subhalo located at dis- 
tance A from the observer along the line of sight defined by 
ip and contributing to the diffuse emission. 
The expression 

K°:r(M,c,r)= f f dcp'dO' f d\' 

J JACl Jl.o.s 



p 2 x (M,c,r(X,\',tP,e'4>')) 
A 2 



J(x,y,z\X' ,0'cp') 



(8) 



describes the emission from such a subhalo. Here, p x (M, c, r) 
is the Dark Matter density profile inside the halo. 

By numerically integrating Eq. 7, we estimate the con- 
tribution to $ cosmo from unresolved clumps in a 1CF 5 sr solid 
angle along the direction ip, for each substructure model con- 
sidered. The result is shown in Fig. 4. In the same figure, 
the solid thick line corresponds to the contribution from the 
MW smooth halo component described in Sec. 2.1 , which is 
computed according to Eq.8, with the distance to the ob- 
server A = Rq . 

We summarize the properties of the smooth subhalo 
contribution in Table 2. In the second and third column we 
show, for each model, the contribution to $ cosmo i n units 
of GeV 2 cm -6 kpc sr towards the Galactic center (ip = 0°) 
and the angle ipd beyond which the smooth subhalo contri- 
bution starts dominating over the MW halo foreground. In 
the fourth column we show the boost factor for each model, 
computed as the ratio of the integral over the MW volume of 
the density squared including subhaloes to the same integral 
for the smooth MW only: 

f M w dV ( P l 



b = 



°MW smooth < rsh) 



- ^ P\lW smooth 



(9) 



To set the two remaining parameters that determine 
the intensity of the annihilation flux, namely the particle 
mass m x and the annihilation cross section cr alm w we adopt 
the most optimistic combination allowed by the SUSY and 
UED models shown in Fig. 3 that do not not exceed the 
current EGRET upper limits for the annihilation flux above 
3 GeV and within a solid angle of 10~ 5 sr. The latter re- 
ceives contribution from two distinct components: the first 
one is of Galactic origin, dominates for ip < 40° and is char- 
acterized by a power-law photon spectrum, that leads, upon 
extrapolation at high energies (Bergstrom et al. 1998) to 
the following parametrization 
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where I and b are the galactic latitude and longitude. The 
normalization factor No depends only on the interstellar 
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Figure 5. Number of photons above 3 GeV, in 1 year in a 
solid angle of 10~ 5 sr, as a function of the angle ip from the 
GC. From top to bottom, the lines correspond to the B Zc 5 CT , 
B Zc , ENSz c , B Z0j5(T , B z „and ENS Z „ model. Empty (filled) circles 
show the B re j Zc (B re y ZQ ) model. ff ann i) = 3 X 10 _26 cm 3 s _1 , 
m x = 40 GeV and BR b ^ = 100% have been used, corresponding 
to the best value <J> PP = 2.6 X 10" 9 cm 4 kpc" 1 GoV^s" 1 sr" 1 . 
The solid thick line shows the EGRET diffuse expected Galactic 
and extragalactic background computed along / = 0. The mass 
function slope is a = 1. 



Results for subhalo models 
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Table 2. Results for the halo models. Column 1: Model 
name. Column 2: <j> cosmo [ GeV 2 cm" 6 kpc sr] toward the 
Galactic Center. Column 3: Angle at which the subhalo 
diffuse contribution dominates over the MW smooth fore- 
ground [degrees]. Column 4: boost factors. Column 5: [ 
10" 9 cm 4 kpc" 1 GcV^s^sr- 1 ]. The value for the zO and the 
B re j Zc models correspond to our best case particle physics sce- 
nario. Values for the zc models except the B re ^ Ze are normalized 
to EGRET data. 
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matter distribution, and is modeled as in Bergstrom et al. 
1998. 

The second one is extragalactic, dominates at ip > 40° , 
and for its photon spectrum we use an extrapolation from 
low energy EGRET data, following Sreekumar et al. 1998: 



7 / extra — y 

"^diffuse = L3 g x 10 -6 £ ,-2.1 



7 



dQ.dE ' cm 2 ssrGeV" ^ 

The thick line in Fig. 5 shows the EGRET photon 
flux as a function of ip. The EGRET flux is computed ac- 
cording to Eqs. 10 and 11 extrapolated above 3 GeV, 
within a field of view of 1CP 5 sr in 1 year of observation. 
The flux is computed along I = 0, that is away from the 
Galactic plane, where the EGRET flux is minimum. The 
Galactic and extragalactic contributions are clearly visi- 
ble. The other curves show the predictions for all mod- 
els in Table 1, obtained when using our best case particle 
physics scenario described in Sec. 3.1 that corresponds to 
$ pp = 2.6 x 1(T 9 cm 4 kpc" 1 GeV" 2 s" 1 sr" 1 . All the values 
of a ann « and ra x that would result in the same <I> PP are 
represented with a solid line in Fig. 3. 

With this choice of parameters all models flagged with 
zc (but the B re / iZc ), for which the halo properties are com- 
puted at the collapse redshift, exceed experimental data 
close to the GC and for ip > 40°, where the EGRET flux is 
assumed to have an extragalactic origin. We do not regard 
the mismatch at small angles as significant, because of the 
limited angular resolution of the original EGRET data. On 
the contrary, we decrease the values of <3> pp to bring the 
zc-models into agreement with data at large angles from the 
GC. The values of $ pp adopted in each model are listed in 
the last column of Table 2. We note that a smaller <J? PP cor- 
responds to assuming a larger particle mass or a smaller 
cross section. The [<7 ann i;, m x ] phase space parameter al- 
lowed for SUSY (circles) or UED (dotted line) models is 
shown in Fig. 3. In the same figure the three dashed and 
dot-dashed lines correspond to the EGRET constrained val- 
ues of $ pp for the zc models. The particle physics models 
above the corresponding line are thus excluded by EGRET 
data. 



3.3 Annihilation flux from individual clumps 

Besides the diffuse signal produced by annihilation in both 
the subhalo population and the smooth MW component, we 
consider here the contribution from individual subhaloes, 
that we regard as Poisson fluctuation of the underlying 
mean distribution of subgalactic haloes that could be de- 
tected as isolated structures. To estimate their flux we con- 
sider 10 independent Monte Carlo realizations of the clos- 
est and brightest subhaloes, in a cone of ~ 50° pointing 
toward the Galactic Center. To do this we generate, for 
each mass decade, the positions of those subhaloes that have 
$ co S mo > j^o^ = 5Q o^ _ 5 x i - 5 GeV 2 cm- 6 kpcsr 

, where the brackets indicate the mean annihilation flux. If 
N < 100 such objects are found, then we still include the 
remaining 100 — N nearest subhaloes in that mass range. 

Adding contribution from an increasing number of indi- 
vidual haloes monotonically increases the chance of subhalo 
detection within the angular resolution element of the de- 
tector. To check whether our procedure is robust and the 
number of detectable haloes has converged we reduced the 



number of Monte Carlo-generated haloes by 70 % and found 
that probability of subhalo detection indeed remains con- 
stant. 

To summarize, for each model, the total contribution to 
jjjcosmo - g gj ven jjy the sum of three terms: the diffuse contri- 
bution coming from unresolved haloes, corresponding to the 
mean contribution of the clumpy component computed with 
Eq.7, the contribution of the diffuse smooth Galactic com- 
ponent and that of individual nearby subhaloes, both com- 
puted using Eq.8. Figs. 6 and 7 show the three contributions 
to $ cosmo for the models B re / iZc and ENS Z0 , respectively. 
In each figure, the contribution from unresolved clumps is 
shown in the upper left panel, the one from the diffuse MW 
in the upper right, the one from resolved clumps in the lower 
left, and the sum of all contributions in the lower left panel. 
The smooth MW halo contribution falls rapidly with the 
distance from the GC, while the diffuse subhalo contribu- 
tion keeps a high value even at large angular distances. The 
single halo contribution is almost completely hidden by the 
overwhelming diffuse foreground, while it is nicely resolved 
as a standalone component. 

The same procedure described in this section for a cone 
pointing to the GC has been repeated for two other regions: 
a cone pointing to the Galactic anticenter and a cone to 
the Galactic pole b = 90°. The purpose is to compute the 
annihilation flux and to evaluate the number of detectable 
subhalos over the whole sky. 



4 PROSPECTS FOR DETECTION 

In absence of strong features in the annihilation spectrum, 
the best chances to detect the annihilation signal within our 
Galaxy is to observe some excess on the 7-ray sky either due 
to diffuse emission or to resolved sources that have no astro- 
physical counterpart. However, the requirements for signal 
detection are different in the two cases since the smooth 
annihilation flux, that contributes to the signal in the first 
case, adds to the noise in the second. 

To determine the probability of halo detection, we con- 
sider a 1 year effective exposure time performed with a 
GLAST-like satellite. We note that, given the 2.4 sr field of 
view and the all-sky survey mode of GLAST, such an expo- 
sure will be achieved in about 5 years of actual observation 
time. 

The prospects for detecting 7-rays from DM annihila- 
tions are evaluated by comparing the number n 7 of expected 
signal photons to the fluctuations of background events n b k g - 
To this purpose we define the sensitivity a as: 



(12) 

/ Af(E, 0;)[#* ignal /dEdQ]dEdQ, 

I Ebk g A lt s ( E > ei)idct> bkg /dEdn]dEdn 



V^bkg 



where Tg defines the effective observation time and 4>bk g is 
the background flux. The quantity £An is the fraction of sig- 
nal events within the optimal solid angle Af2 corresponding 
to the angular resolution of the instrument and is usually 
~ 0.7. We set it equal to 1 to get the most optimistic val- 
ues. The effective detection area A eB for electromagnetic or 
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Figure 6. Map of <J> cosmo (proportional to the annihilation signal) for the B re j- 2c model, in a cone of 50° around the Galactic Center, as 
seen from the position of the Sun. Upper left: smooth subhalo contribution from unresolved haloes. Upper right: MW smooth contribution. 
Lower left: contribution from resolved haloes. Lower right: sum of the three contributions. 



hadronic particles is defined as the detection efficiency times 
the geometrical detection area. In the following we make 
the realistic assumption that all hadronic particles will be 
identified, so that the background will be composed by pho- 
tons only. We also assume that all photons will be correctly 
identified, which is somehow optimistic, since there will be 
a small amount (a few percent) of irreducible background 
due, e.g., to the backsplash of high energy photons. 
We use A bS = 10 4 cm 2 , independent from the energy E and 
the incidence angle $i, and an angular resolution of 0.1°. 
Both these values are rather optimistic, since the expected 
GLAST angular resolution approaches 0.1° only at about 20 
GeV, while the on-axis effective area is quoted to be max- 
imum ~ 8 x 10 3 cm 2 above 1 GeV and decreases by ~20% 
for an incidence angle of 20°. 

As anticipated, different annihilation signals need to be 
compared with different background noises. For the detec- 
tion of the diffuse annihilation flux the background is con- 
tributed both by Galactic (Eq. 10) and extragalactic astro- 
physical sources (Eq. 11) measured by EGRET. In the case 
of individual subhaloes, the annihilation photons produced 
in the smooth Galactic halo and in the unresolved clumpy 
component contribute to the background rather than to the 
signal. 



4.1 Sensitivity to diffuse emission 

We first study the sensitivity a of such a GLAST-like obser- 
vatory to the annihilation flux from the smooth DM profile 
and from the diffuse contribution of unresolved subhaloes. 
Both signals are computed above 3 GeV, and the astrophys- 
ical background noise is obtained from Eqs.10 and 11 spec- 
ified along / = 0. 

The result is shown in Fig. 8, where we plot the statisti- 
cal significance of the detection as a function of ip for each of 
the models listed in Table 2.2. 1 a detections of the annihila- 
tion signal is expected at ip < 40° for all the models labelled 
zc. The chances of observing the diffuse annihilation flux are 
significantly higher in the direction of the Galactic Center 
along which models labelled zO predict a signal detectabil- 
ity as high as 5 a. Yet, these predictions should be taken 
with much care since the measured astrophysical 7-ray flux 
above 3 GeV in the direction of the GC, which constitutes 
the background, is known with large uncertainties. 

4.2 Detection of individual haloes 

Subhaloes can also be detected through the annihilation flux 
produced by individual, nearby clumps that would appear as 
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Figure 7. Same as in Fig. 6 for the ENSz model. 



bright, possibly extended, sources, as shown in the bottom 
right panels of Figs. 6 and 7. In this case the signal is pro- 
duced within the individual haloes of our Monte Carlo real- 
izations, while the background is contributed by the smooth 
astrophysical background plus the diffuse annihilation flux 
produced by the Galaxy and its subhaloes. 
For each halo in the 10 Monte Carlo realizations, and for 
each virial concentration model, we 

• assign to it an arbitrary concentration parameter c(M); 

• calculate the annihilation signal; 

• find the value of the concentration parameter that guar- 
antees a 5 a detection in 1 year exposure time, C5 CT (Af); 

• identify the probability of detection of the clump with 
the probability P(> cs a ) that such a clump has a concentra- 
tion as high as cs a (M), assuming the lognormal distribution 
described in Eq.4. 

The total number of detectable subhaloes is then sim- 
ply given by ^\ Pj(> cs„), where the sum is performed over 
all haloes in the realization. Results are obtained by aver- 
aging over all over the 10 Monte Carlo realizations and the 
procedure is repeated for all models listed in Table 1. 

The number of haloes that can be detected in 1 year 
with a significance above 5 a in cone of view with angular 
opening of 50° towards the GC is shown in Fig. 9 for the zO 



models and in Fig. 10 for the zc models, as a function of the 
subhaloes mass. Had we assumed a deterministic relation 
for c(M) the number of events would have decreased by a 
factor ~ 2. 

At higher latitudes the number of haloes indeed reduces, 
but this is compensated by the lower foreground given by 
the smooth subhalo component. The maximum number of 
detectable events is obtained toward b = 90 where these two 
effects interplay in a most favourable way for the detection. 

Table 3 lists the number of haloes that can be detected 
with a significance larger than 5 a in a cone of 50° around 
the Galactic Center (first column), the Galactic pole (sec- 
ond column) and the Galactic anticenter (third column), for 
each model, for our reference mass function slope a = 1 and 
the $pp values listed in Table 2. Hereafter, each error is the 
standard deviation obtained averaging over the 10 MC rep- 
resentations. The effect of decreasing the number of haloes 
far from the GC is compensated by the lower foreground due 
to the diffuse subhalo contribution to the annihilation flux. 
The best compromise is found around the Galactic poles. 

In Table 4 we show the total number of haloes that can 
be detected with a GLAST-like satellite in the whole sky 
with a significance larger than 5 a, with a mass function 
slope a — 1 (first column) and &pp values listed in Table 2. 
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Figure 8. Statistical significance, as a function of the angle 
from the GC tp, for the detection of the DM annihilation flux 
from diffuse subhaloes plus the MW smooth component (along 
1 = 0), for the different models explored: B ZQ (dotted line ), 
ENS ZQ (short dashed), B re ^ Z[| (filled circles), B re ^ Zc (empty cir- 
cles), B ZOj 5 CT (long dashed), B Zc (long dot-dashed), ENS Zc (dot- 
dashed), B Z( . i 5 (T (solid). The mass function slope is a = 1. We 
refer to Table 2 for the values of <J>pp used in this figure. 



Number of detectable haloes (a = 1) 



Model 


ly GC 


m5<t 
JV 90 


JV5ct 
IV 180 


B r e/, Z0 


0.65 ±0.45 


0.85 ± 0.43 


0.59 ±0.30 


B Zf) 


0.65 ±0.45 


0.84 ± 0.43 


0.59 ±0.30 


B Z0 ,5ct 


0.46 ±0.34 


0.60 ±0.34 


0.50 ±0.27 


B r e/, Zc 


16.16 ±2.60 


23.24 ±2.28 


18.55 ± 1.72 




1.74 ±0.92 


2.40 ±0.82 


2.15 ±0.65 


B Zc ,5ct 


0.05 ± 0.05 


0.07 ±0.08 


0.08 ± 0.08 


ENS Z0 


0.06 ±0.12 


0.07 ±0.10 


0.04 ±0.07 


ENS Zc 


0.29 ± 0.40 


0.49 ± 0.37 


0.46 ± 0.30 



Table 3. Number of haloes detectable, at 5 a in 1 year of effec- 
tive observation with a GLAST-like satellite, in a 50° f.o.v. cone 
towards the GC (column 1), the Galactic pole (column 2) and 
the anticenter (column 3). The subhaloes mass function slope is 
a = 1. We refer to Table 2 for the values of Qpp used in this 
table. 



The remaining two columns show the effect of adopt- 
ing a mass function with power-law index a = 0.95 (sec- 
ond column) and a = 0.9 (third column). Adopting a shal- 
lower mass function increases in most cases the number of 
detectable subhalos. However, the magnitude of the effect, 
that results from the lowering of the unresolved background, 



Total number of detectable haloes 
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B t 

*-*rej , Zc 


118.36 ±24.96 
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125.03 ± 20.06 


B, c 
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104.23 ± 24.78 


119.04 ±19.77 


B Zc ,5ct 


0.39 ±0.56 


10.55 ±6.36 


96.34 ±18.66 


ENS Z0 


0.33 ±0.89 


0.67 ± 1.58 


0.34 ±0.50 


ENS Zc 


2.50 ± 4.48 


23.43 ± 10.17 


30.40 ±10.31 



Table 4. Total number of haloes detectable over the whole sky, at 
5 <r in 1 year of effective observation with a GLAST-like satellite, 
for a mass function slope a = 1 (column 1), a = 0.95 (column 
2) and a = 0.9 (column 3). We refer to Table 2 for the values of 
<I>PP used in this table. 

whose main contribution is given by small sub-haloes, de- 
pends on the model explored. 

As expected, we can observe how this effect is larger 
for those models whose overall contribution to <J> cosmo ; s 
larger, that is for those models whose concentration param- 
eters have been computed at the collapse redshift. The effect 
is reduced for the other models, as well as for the B re / jZc one, 
for the following reason: when using a = 1, the 3> pp value 
for the zc models (but the B re / iZe one) has been decreased 
in order to respect the EGRET EGB limit, while when us- 
ing a = 0.9 all models fulfill the EGRET EGB constraint, 
and we can use our best case $pp. This is also true for 
q = 0.95, except for the B Zc : 5 CT model, where we have to use 
$ pp = 0.84cm 4 kpc^GeV^s^sr- 1 . In fact, the zO and 
the B re / jZc model experience just a minor increase (compati- 
ble within the error bars) of the number of detectable halos; 
the B Z( .and ENS Zc models reach this stability for a < 0.95, 
as soon as their $ pp allowed value gets our best value; the 
B z< .,5ct model keeps on showing a large effect when changing 
mass function slope, because it is allowed to have the best 
value <I> PP only when a = 0.9. 

In the most optimistic B re / jZc model, we expect that a 
GLAST-like experiment could detect ~ 120 — 130 subhaloes 
with masses above 10 5 M© over all sky, for all the mass 
function slopes considered in this analysis. In all the models 
whose concentration parameters are computed at z = the 
number of detectable events is compatible with zero within 
the errors, whatever slope is used. Accordingly to the afore- 
mentioned discussion, the effect of changing the mass func- 
tion slope is dramatic in the B Zc (B Zei 5 CT ) model, for which 
the total number of events ranges from ~ 10 (~ 0) for a — 1 
to ~ 120 (~ 100) for a = 0.9. A large effect (~ to ~ 30) 
is observed in the ENS Zc model too. 

5 DISCUSSION AND CONCLUSIONS 

The prospects for detecting 7-rays from the annihilation of 
DM particles in substructures of the MW have been inves- 
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tigated by a number of authors (e.g. Stoehr et al. 2003, 
Pieri & Branchini 2004, Koushiappas et al. 2004, Oda 
ct al. 2005, Pieri et al. 2005). In this work we confirm 
that substructures can provide a significant contribution 
to the expected Galactic annihilation signal, although the 
actual enhancement depends on the assumptions made on 
the clump properties, which are affected by large uncertain- 
ties. Indeed, given the assumed substructure mass function 
dN/dM oc M~ 2 , the contribution to the total 7-ray flux 
by subhaloes of different masses depends on the annihila- 
tion signal produced within each clump which is dictated by 
the internal structure. Numerical experiments have shown 
that the total annihilation signal is dominated by the high- 
est mass subhaloes both in a Galactic halo at z=0 (Stoehr 
et al. 2003) and in 0.1 Mq halo host at z=75 (Diemand et 
al. 2006). However, the recent results of the high resolution 
'Via Lactea' simulation (Diemand et al 2007a) indicate that 
the annihilation luminosity is approximately constant per 
decade of substructure mass while analytical calculations 
(Colafrancesco et al. 2006) tend to find that the signal is 
dominated by small mass subhaloes. Indeed this is also the 
case with our model predictions. Fig. 11 shows the expected 
contribution of the unresolved haloes to the total annihila- 
tion flux as a function of the subhalo mass, integrated on 
each mass decade. In all the models explored the annihila- 
tion signal is dominated by the smallest clumps, as a result 
of the decrease of the virial concentration with the subhalo 
mass, as shown in Fig. 2. Under optmistic assumptions on 
the particle physics parameters of DM particles, a GLAST- 
like experiment might detect such a DM annihilation flux, 
but only in a few pixels around the Galactic Center. 

It should be noticed, however, that estimates of the an- 
nihilation signal from the Galactic Center are affected by 
the poor knowledge of the DM profile in the innermost re- 
gions of the Galaxy, which is usually obtained by extrap- 
olating over many orders of magnitude the results of nu- 
merical simulations. The presence of a Supermassive Black 
Hole at the center of the Galaxy makes things even more 
complicated, as it may significantly affect the distribution 
of DM within its radius of gravitational influence , leading 
to the formation of an overdensity called "spike" (Gondolo 
and Silk, 1999). Spikes require however rather fine-tuned 
conditions to form (Ullio et al. 2002) and any overdensity is 
in any case severely suppressed by the interaction with stars 
and DM self-annihilations (Merritt et al. 2002, Bertone & 
Merritt 2005). 

In alternative, one could look for an annihilation signal 
from individual DM substructures, such as dwarf galaxies 
or even smaller, 'baryon-less', clumps. We have shown that, 
depending on the assumptions made on the properties of 
clumps, only large haloes with M > 1O 5 M0 can be detected 
and identified with a GLAST-like experiment, which is con- 
sistent with the analyses of Stoher et al. 2003 and Koushi- 
appas et al. 2004. The number of detectable haloes ranges 
from to more than a hundred, depending on the model. 
Adopting a shallower subhalo mass function increases the 
number of detectable subhaloes in those models in which the 
diffuse annihilation signal is dominated by the unresolved, 
low mass haloes. 

In any case, scenarios leading to a large number of 
detectable small-scale clumps appear to be severely con- 
strained by the 7-ray background measured by EGRET. 
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Figure 9. Number of events detectable in 1 year a 5 a with a 
GLAST-like experiment in a 50 degrees cone towards the GC for 
the models ENS Z0 (dot-dashed), B ZQ (solid), B Z0 .5 CT (dashed) and 
B r e/,z (l° n S dot-dashed), assuming the lognormal distribution 
for the concentration parameter. We refer to Table 2 for the used 
values of $pp. The mass function slope is a = 1. 
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Figure 10. Number of events detectable in 1 year a 5 ct with a 
GLAST-like experiment in a 50 degrees cone towards the GC for 
the models B re ^ jZe (dashed) ENS Zc (solid), B Zc (long dot-dashed) 
and B Z( . j 5 (T (dot-dashed), assuming the lognormal distribution for 
the concentration parameter. We refer to Table 2 for the used 
values of $pp. The mass function slope is a = 1. 
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In particular, the model of Koushiappas 2006 is similar 
to our B Zc model, as far the cosmological term is considered, 
while the particle physics contribution corresponds to our 
best case scenario <&gf . Although nearby haloes would be 
bright, and observable, in this case, we have shown that the 
associated diffuse emission produced by all the remaining, 
unresolved, clumps in the Milky Way, would far exceed the 
7-ray background measured by EGRET. The chances to de- 
tect the proper motion of clumps are thus very low, as the 
lowest mass detectable subhaloes, (M=10 Mq), are typi- 
cally found at a distance greater than 0.5 kpc, leading to a 
proper motion less than ~ 0.1' yr _1 , well below the GLAST 
angular resolution of a few arcminutes. 

We have made use of simplified and extreme scenar- 
ios for the subhaloes concentration parameter models. More 
accurate scenarios, though not supported by numerical sim- 
ulations for small mass haloes, could lead to different dif- 
fuse foreground levels and to both more or fewer detectable 
haloes. 

One may wonder why we preferentially expect to indi- 
vidually detect the more massive subhaloes, while the un- 
resolved annihilation signal is mainly contributed by small 
mass clumps. The reason is that the volume over which 
individual haloes can be detected decreases rapidly with 
the halo mass. To see this, let us consider the maximum 
distance Umax at which a clump can be detected. This 
distance depends on the halo luminosity which, in turns, 
depends on the halo mass and concentration. For a NFW 
profile: Dmax oc M°' 5 c(M) 15 . On the other hand, as dis- 
cussed e.g. by Koushiappas 2006, given a subhalo mass func- 
tion dN/dM oc M -2 , the number of detectable haloes per 
mass decade is: dN/dLog(M) oc D^^M' 1 oc M°' 5 c(M) 4 - 5 . 
Assuming a simple scale-free virial concentration c(M) oc 
M ', we see that if 7 > 7^ = —1/9, then the number 
of detectable haloes per mass decade halo indeed increases 
with the subhalo mass. The c(M) relation for some of our 
models is shown in Fig. 2 along with the 7^ = — 1/9 refer- 
ence slope (thick dashed line). Models B zo and ENS Z „, that 
have 7 > jth do indeed predict that the probability of sub- 
halo detection increases with the mass (see Fig. 9). On the 
other hand, the slope of the c(M) relation for M < 1O 2 M0 
is sligthly steeper than 7^. Therefore we expect a bimodal 
probability that peaks at high masses with a secondary max- 
imum for the smallest subhaloes. Indeed, this is what we 
observe in Figure 9. 

In conclusion, we have studied the prospects for indirect 
detection of Dark Matter in MW subhaloes with a GLAST- 
like satellite. We have chosen 8 different models for the con- 
centration parameter, which span the phase space the the- 
oretical uncertainties on the Dark Matter halo properties, 
as well as 3 different values for the subhaloes mass function 
slope. For each model, we have computed the diffuse emis- 
sion from unresolved subhaloes, as well as the 7-ray flux 
from individual, nearby haloes. 

We found that for models with concentration param- 
eter computed at z — the detection of individual haloes 
appears challenging, while the diffuse emission from unre- 
solved clumps dominates the MW smooth emission for sky 
directions > 1 degree off the GC. 

In the case of a = 1, in all the zc models except the 
B re / jZe , the diffuse emission from unresolved clumps exceeds 
the EGRET constraints in a portion of the DM parame- 
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Figure 11. Contribution to Qcosmo from substructures of mass 
M integrated over the mass decade, for the different B mod- 
els explored and computed toward the Galactic Center. Solid 
and empty circles correspond to the B re j Z[| and B re ^ Zc models 
respectively. Lines show the B Zci 5 CT (dot-dashed), B Zc (dashed), 
B Z0 ,5 CT (dotted) and B zo (solid) models. The mass function slope 
is a = 1. 

ter space relevant for SUSY models, as shown in Fig. 3. 
Adopting DM models compatible with the EGRET data, 
one may still hope to detect individual haloes, like e.g. in 
the B Zc model. 

The B re / iZc is our best case model for all the values of 
the mass function slope, though it should be stressed that 
it is not supported by numerical experiments but only by 
theoretical considerations. The success of the B re /. Z( . model 
is due to the fact that the diffuse emission expected from 
subhaloes is dominated by small mass haloes while the large 
mass haloes are most favourably detected as spare sources. 
A functional form for the concentration parameter which 
flattens at low masses, as it is the B re / jZc one, will in fact 
decrease the diffuse emission, thus allowing a larger value for 
<I>pp and consequently increasing the chances of detection 
for large mass haloes. 

In general, adopting the most optimistic set of parame- 
ters for the DM particle compatible with EGRET (see Table 
2), the number of detectable subhaloes over all sky (at 5a, 
with a GLAST- like experiment and a 1-year exposure time), 
for the 8x3 models we have studied, ranges between and 
120. 

Yet, it should be noticed that the numbers listed in 
Table 2 are obtained with a very optimistic value for the 
Particle Physics involved in the process. If we assume a 
more realistic model for the Dark Matter particle ( m x = 
100 GeV , av = 10 -26 cm 3 s _1 ) instead, we find that at most 
only a handful of detectable haloes are found (at most a 
handful over all sky for the most optimistic B re / jZc model). 

In all the models explored, small mass subhaloes are al- 
ways below the threshold for detection, and their presence 
could be revealed only through the enhancement of the dif- 
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fuse foreground emission. The different predicted ratio be- 
tween diffuse emission and number of detected haloes for 
the models we have considered, would provide precious in- 
formation on the underlying cosmology, in case of positive 
detection. 
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